Human health risk assessment of heavy metals in vegetables of Bangladesh

This study aims to evaluate the heavy metal concentration in fifteen species of vegetables as well as associated health risk. Atomic absorption spectrometry is used to assess heavy metals. The mean concentrations of Pb, Cd, Cr, Ni and Fe in vegetables were 4.78, 0.713, 9.266, 0.083, 5.06 mg/kg/fw exceeding the reference value of FAO/WHO indicating unsafe to consumption. Based on principal component analysis, the Pb, Cr, Ni and Fe are from same sources. Health risk was estimated in terms of estimated daily intake (EDI), target hazard quotient, hazard index (HI) and cancer risk (CR). The EDI values of metals except Cr were found to be lower than maximum tolerable daily intake (MTDI). The total THQs of metals were > 1 indicating non-carcinogenic health risk. The individual HI values for vegetables except potato (0.831) and total HI values were found to be > 1 (94.747). The TCR of Pb, Cd and Cr were > 1.0E−04 which indicating carcinogenic risk. Fruit and pod vegetables contribute much in carcinogenic risk for Pb and Cr whereas fruit, root and stems vegetables for Cd. The study revealed potential human health risk associated with the consumption of different types of vegetables in Bangladeshi adult population that might assist the regulatory bodies to develop new strategies to minimize the risk to human.

www.nature.com/scientificreports/Noakhali is a coastal region of Bangladesh where different types of industries are established.Heavy metal pollution in coastal region is very common due to discharge of heavy metals in coastal water by several anthropogenic activities such as food processing, pharmaceuticals, paper industry etc. 16,17 .This have significant consequences on environment due to bio-magnification in food chains 18 .However, there is no specific data of heavy metals contamination in vegetables grown in the area of Noakhali.As there is lack of information about heavy metal contents in most common consumed vegetables in Noakhali, Bangladesh as well as differences of metal concentration among vegetable types, the study aims to represent metals concentration in vegetables and potential carcinogenic and non-carcinogenic health risk for these vegetables' consumption.The study also evaluates the sources and differences of heavy metals in the vegetables.

Methodology Sampling and Study area
The study was carried out in Noakhali district (N22°49′28.7″,E91°6′6.24″) of Bangladesh during the year of 2022.Fifteen different sample species of vegetables were collected from three different local markets located in Noakhali district.The sampling sites were shown in Fig. 1 which was created by QGIS software (version 3.10.2).For each species, 3 replicates vegetable sample were collected from each market.In total, 135 samples were collected from three local markets.All samples were washed and stored in fresh polybag and brought to laboratory for analysis.

Wet digestion
In the sample preparation, all samples were cleaned with deionized water and chopped with stainless steel blade.Then, all samples were dried in oven at 100 0 C to remove all moisture and grinded with the help of a mortar.Wet digestion method was used with a slight modification of the method by Akinyele et al. 19 .To digest the samples, acid mixture was used (69% concentrated HNO 3 :95-97% concentrated H 2 SO 4 = 1:4).According to digestion, 0.5 g of each homogenous dried sample was placed in digestion tube and added 5 ml of acid mixture into it.Then the mixture was heated from 130-170 °C for an hour until transparent solution was observed.After cooling, 2 ml of H 2 O 2 (30% concentrated) was added and heated until clear solution was found.Then the clear solution was filtered using Whatman 102 filter paper and then diluted to 250 ml of deionized water.Then, the samples were collected in falcon tube and stored for analysis by atomic absorption spectrometry.

Atomic absorption spectrometry (AAS)
The heavy metal determination was carried out using a Perkin-Elmer novAA 900H atomic absorption spectrometer 20 .To detect each metal, a single beam hollow cathode lamp of chromium, cadmium, lead, iron,

Quality control
To ensure the quality control, blank sample was used in every sample analysis.The reliability of each method was evaluated by linearity, range and recovery.The procedure of validation of result were summarized in Supplementary Tables 1 and 2. The percentage of recovery (% R) was ranged from 97.15 to 101.55% and the linearity was ranged from 0.995 to 0.998.Effective aseptic technique was followed to prevent contamination in samples and all glass materials were rinsed with 1% nitric acid.

Estimated daily intake (EDI) of heavy metals
EDI was measured in mg/kg-weight by following formulae 21 .
Where MC is the metal concentration in vegetables (mg/kg fw), IR (g/day/person) is the ingestion rate of vegetables adopted from household income and expenditure survey (HIES) 13 for an adult individual of 60 kg (adult) body weight 22 .

Non-carcinogenic risk
The non-carcinogenic risk of heavy metals due to consuming vegetables were estimated in terms of target hazard quotient (THQ), total target hazard quotient (TTHQ), and hazard index (HI).THQ THQ was calculated as per USEPA Region III Risk-based Concentration Table 23 and Wang et al. 24 .The equation used for estimating THQ was; where EF is the exposure frequency (365 days/year); ED is the exposure duration (70 years for non-cancer risk in this study), as used by 25 ; FIR is the food ingestion rate (g/person/day); CM is the heavy metal concentration in vegetable (mg/kg); BW is the average body weight (bw) (adult: 60 kg); AT is the average exposure time for non-carcinogens (EF × ED) (365 days/year for 70 years; i.e.AT = 25,550 days); RfD is the oral reference dose of the metal.RfD values of Cd, Cr, Ni, Pb and Fe are 0.001, 1.5, 0.02, 0.0035 and 0.667 (mg/kg/day), respectively 26 .
If the THQ < 1, the exposed consumers are unlikely to experience any adverse health risk, while if the THQ ≥ 1, there is a potential health risk 24 , and associated interventions and protecting initiatives are required to be taken.TTHQ TTHQ for individual from THQs is expressed as the sum of the hazard quotients 25 .
where TTHQ < 1 is safe, TTHQ > 1 is hazardous, and THQ (Cr) is the target hazard quotient for Cr intake.

Hazard risk (HI)
HI is assessed to estimate overall potential for non-carcinogenic health risk from consuming more than one metal.

Target cancer risk
The method to estimate TCR is also provided in USEPA Region III Risk-Based Concentration Table 25 .The model for estimating TCR was shown as follows (4): where EDI is the estimated daily intake, CPSo is the carcinogenic potency slope for oral route of 0.0085 (mg/ kg-day) −1 for Pb, 15 (mg/kg-day) −1 for Cd, and 0.42(mg/kg-day) −1 for Cr 27 .In general, CR value lower than 1.0E−06 is considered to be negligible, above 1.0E−04 is considered unacceptable, and lying between 1.0E−06 and 1.0E−04 is considered an acceptable range 28,29 .Data required for estimating EDI, THQ and CR are summarized in Supplementary Table 3.

Ethical consideration
The study was classified as exempt according to the institutional ethics committee of the Noakhali Science and Technology University.All methods were performed in accordance with the relevant guidelines and regulations.

Statistical analysis
All statistical analyses were performed with SPSS 23.0 Inc., Chicago, IL, USA for Windows.Data were presented as mean and standard deviation (SD) and were subjected to one-way analysis of variance (ANOVA) to assess whether heavy metals varied significantly between vegetables (p < 0.05).To check the similarities and differences of heavy metals distribution, principal component analysis (PCA) and cluster analysis (CA) were performed.
The PCA analysis was done using Varimax normalized rotation method to maximize the sum of the variance of the factor coefficient and CA was performed using Ward's method.

Concentration of heavy metals in commonly consumed vegetables
Heavy metal concentration (Pb, Cd, Cr, Fe and Ni) in vegetables (mg/kg fw) were presented in Table 1.The concentration of heavy metals was increasing in the following order of Fe < Cd < Pb < Ni < Cr.The mean concentration of metals in vegetables was decreasing in order of pointed gourd > yard-long bean > cauliflower > teasel gourd > egg plant > ladies finger > taro > bitter gourd > Chinese okra > snake gourd > carrot > pumpkin > tomato > papaya > potato.Highest concentration of Cr was observed in vegetables and in Bangladesh, the sources of Cr are utilization of industrial and untreated water in vegetable production, use of chemical fertilizer and pesticides 14,30 .Compared with other vegetable types, the concentrations of Pb, Cr, Fe and Ni were higher in pod vegetables.Although studies found that leafy vegetables have high accumulation of heavy metals as they absorb large amount of heavy metals from soil and have large surface areas 31 .The concentration of heavy metals in different types of vegetables depends on the soil properties and accumulation capacities of vegetables 32 .The concentration of metals in vegetables were higher than the recommended value set by FAO/WHO except Fe.
In the case of Fe, it is a naturally occurring element that is essential for plant growth and development.Plants have evolved mechanisms to regulate Fe uptake and distribution within their tissues to meet their physiological requirements while avoiding toxicity 33 .The highest concentration of Pb was found in cauliflower and yardlong bean.The mean concentration of Pb, Cd, Cr and Ni were higher than the previously studied literature in Bangladesh and Pakistan 4,7,34-36 (Table 2) whereas the mean concentration of Fe in vegetables was higher than Ahmed et al. 34,36 and similar to Sultana et al. 37 (Table 2).Iron is an essential metal to prevent anemia however, excess consumption of iron may cause heart disease and liver damage 38 .

Principal component analysis and cluster analysis
Varimax-normalized rotation method was used for principal component analysis for estimating the factor loadings in each metal (Fig. 2).According to the results, two eigen values were greater than one and the first two component were explained the variance by 39.63% and 26.21% (Supplementary Table 4).PC1 revealed highest loadings for Pb, Cr, Ni and Fe which indicated that their source of origin was same and mostly contributed by anthropogenic activities such as use of fertilizers, pesticides, organic matters etc. 39,40 .Further, we used cluster analysis (CA) using Ward's method with dendrogram for dividing the vegetables into different species (Fig. 3).Mean concentration of heavy metals was used for CA.Different cluster was formed between the vegetable species according to their similarity of nature.www.nature.com/scientificreports/

Dietary intake of heavy metals
Assessment of dietary intake of food is an essential tool for measuring the amount of nutrients intake which may be lead to deficiencies or health risks 41 .According to the average concentration of heavy metals in vegetable, dietary intake of metals for adult individual was estimated and presented in Table 3.Total dietary intake of metal of Pb, Cd, Cr, Fe and Ni from vegetables were 0.185, 0.0255, 0.3823, 3.4059, and 0.2057 mg/day for adults respectively.Dietary consumption of all of the metals except Cr from vegetables were less than the value of maximum tolerable dietary intake (MTDI) set by FAO/WHO which indicates that there is no potential human risk prior to consumption of vegetables from the studied area.This study showed that the contributions of studied vegetables to daily intake of Pb, Cd, Cr, Fe and Ni were 88%, 55%, 190%, 20% and 70% of the maximum tolerable dietary intake.

Non-carcinogenic and carcinogenic health risk
The non-carcinogenic risk in terms of target hazard quotient (THQ), total target hazard quotient (TTHQ) and hazard index (HI) and carcinogenic risk of consuming studied vegetables were presented in Table 4.The THQ values of Pb and Cd for almost all species of vegetable were higher than 1, indicated that consumption of these vegetables might cause non-carcinogenic risk.The THQ values of the studied vegetables (except Cr) were higher than 1 which indicated that people might have potential health risk prior to vegetables consumption.A study in India revealed potential health effects from Cd and Ni due to the value of HQ exceeded the safe limit 42 .The ranking order of total THQ for vegetables species were cauli flower > taro > pointed gourd > yard-long bean > eggplant > Chinese okra > snake gourd > teasel gourd > pumpkin > ladies finger > bitter gourd > tomato > carrot > papaya > potato.Considering all the metals, the total THQ (sum of individual metals or HI) was 35.36, 32.44, and 27.017 for fruit vegetable, pod vegetable and root and stem vegetable respectively, in a total 94.797 which is > 1 indicating that these vegetables were not safe for human consumption and consumption on regular basis is not recommended.It is stated that exposure to two or more heavy metals can cause additive or interactive effects 24 .
As Pb, Cd and Cr contributed both non-carcinogenic and carcinogenic health risk depending on the exposure dose and duration, the present study estimated the TCR values of Pb, Cd and Cr due to exposure from several vegetables presented in Table 4. TCR values for Pb ranged from 3.8E-06 to 2.0E-04; for Cd, TCR ranged from 2.3E−03 to 8.3E−02 whereas 1.95E−03 to 1.93E−02 for Cr in vegetables.The TCR values of Pb for pointed gourd, yard-long bean and cauliflower were higher than 1.0E−04 29 which indicated that potential carcinogenic risk of the consumers in the studied area.TCR values of Pb for other vegetable except potato were lying between the range of 1.0E-06 to 1.0E-04 which are considered as acceptable 29 .The TCR of Cd and Cr from the consumption of vegetables were higher than USEPA acceptable limit risk (10 -4 ) indicating that consumers of these vegetables www.nature.com/scientificreports/are exposed to Cd and Cr with a life-time cancer risk.The percentage of carcinogenic risk of Cr (49.61%) and Cd (42.22%) were higher among consumers who intake fruit vegetables whereas higher carcinogenic risk for Pb was found in pod vegetables (42.93%) consumption (Fig. 4).The cumulative cancer risk percentage of Pb (24.6%) and Cr (16.66%) from roots and stems consumption were lower.
In the present study, we have several limitations.We didn't measure the accumulation rate of each vegetable from soil to edible parts.We didn't account the other sources of heavy metal exposure such as dermal, nasal which may associate with health risk to humans.The effect of environment such as seasonal variations on heavy metals concentration in vegetables also need to be addressed.Given the well-established presence of arsenic contamination in Bangladesh 43 and its potential impact on agricultural produce 43 , the omission of arsenic analysis in this study is a significant shortcoming.Limited resources and analytical capabilities for arsenic detection contributed to the exclusion of arsenic analysis in this study as arsenic analysis in vegetables typically requires hydride generation system and the additional sample preparation steps involved in arsenic analysis compared to the detection of heavy metals like lead, cadmium, chromium, nickel, and iron 43 .Another limitation is that the AAS technique used in this study measures the total chromium content, which includes both Cr (III) and Cr (VI), without distinguishing between the different oxidation states.In the absence of speciation data, applying the CPSo value for Cr (VI) to total chromium concentrations may overestimate the carcinogenic risk, as it assumes that all the chromium present in the samples is in the hexavalent form.

Conclusion
The present study revealed the concentrations of heavy metals in commonly consumed vegetables in Noakhali district of Bangladesh and assessed the potential health risk prior to vegetable consumption in terms of THQ and TCR.In Noakhali, maximum of vegetables had toxic metals (Pb, Cd, Cr and Ni) which are higher than the maximum allowable concentration (MAC) however, in case of EDI, only Cr exceed the permissible limit (MTDI) set by FAO/WHO.The total THQ for Pb, Cd, Fe and Ni were > 1 through vegetables consumption indicating potential health risk.Consumption of fruit vegetables, pod vegetables and root and stem vegetables would be unsafe as HI value was > 1. Considering TCR values, the total TCR values for Pb, Cd and Cr were higher than 1.0E−04 suggesting potential health risk from vegetable consumption.The findings emphasize the need for immediate action to mitigate these risks and ensure food safety.Future research should focus on identifying the sources of contamination, assessing other exposure routes, and exploring the influence of environmental factors.Policymakers and stakeholders must collaborate to develop and implement stringent regulations, monitoring programs, and techniques to address this pressing issue.Public awareness campaigns are crucial to educate farmers and consumers about the risks and promote safe agricultural practices.

Figure 2 .
Figure 2. Principal component analysis of heavy metals in vegetables by varimax normalized rotation method showing loading of metals.

Figure 3 .
Figure 3. Cluster analysis (CA) of the vegetables from Noakhali region of Bangladesh.

Table 1 .
Concentration of heavy metals (mg/kg fw) in commonly consumed vegetables in Bangladesh.

Table 2 .
Comparison of average heavy metals (mg/kg fw) in vegetables with previous studies in Bangladesh.

Table 3 .
Comparison of estimated daily intake (EDI) of heavy metals from vegetable with maximum tolerable daily intake (MTDI).

Table 4 .
Carcinogenic and non-carcinogenic risk of heavy metals consuming vegetables in Bangladesh.